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ABSTRACT 


An  analysis  and  test  rig  measurements  of  the  per 
formance  of  a transonic  axial  turbine  are  reported.  The 
purpose  was  to  confirm  the  accuracy  of  measurements  made 
in  a test  rig  which  was  designed  to  separate  the  losses 
occurring  in  the  stator  from  the  losses  occurring  in  the 
rotor  blade  rows.  The  analysis  was  programmed  for  the 
Hewlett-Packard  21 -MX  computer.  Reasonable  agreement 
between  predicted  and  measured  characteristics  was  obtained 
using  experimentally  determined  losses  in  the  computer 
program.  Lack  of  agreement  was  noted  using  theoretical 
values.  It  was  concluded  that  the  rotor  was  not  choked 
at  the  conditions  in  the  tests,  and  that  the  test  rig 
measurements  were  valid.  A successful  technique  for 
smoothing  the  data  obtained  from  the  rig  is  also  reported. 
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I.  INTRODUCTION 


The  transonic  turbine  test  rig  installation  at  the 
Naval  Postgraduate  School  Turbopropulsion  Laboratory  was 
designed  to  study  the  effects  on  turbine  performance  of 
varying  axial  and  tip  clearances,  to  study  the  effects  of 
blading  design  on  turbine  performance,  and  to  allow  the 
separate  determination  of  stator  and  rotor  losses  in  an 
operating  machine. 

Until  the  work  of  Solms  (Ref.  1)  the  separation  of 
rotor  and  stator  losses  through  test  rig  measurements  had 
not  been  attained  satisfactorily.  Through  improvements 
in  hardware  and  instrumentation  and  improvements  in  the  data 
reduction  process  the  stator  and  rotor  losses  were  deter- 
mined separately  and  were  reported  in  Ref.  1. 

Anomalies  remained,  however.  Specifically,  the  turbine 
configuration  designated  as  Turbine  C in  Ref.  1 gave 
different  results  when  compared  in  terms  of  "referred” 
quantities  depending  on  whether  the  discharge  was  to 
atmospheric  pressure  or  to  a region  of  reduced  pressure. 
(Turbine  C had  converging-diverging  stator  passages  in  an 
axial  entry,  single  impulse  stage.  The  turbine  was  designed 
to  operate  in  the  transonic  range.)  In  addition,  R&f.  1 
reported  considerable  scatter  in  the  loss  coefficients. 

In  the  work  of  Robbins  (Ref.  2)  it  was  shown  that 
discharge  pressure  affected  the  measurements  of  flow  rate 
into  the  stage  and  also  affected  the  labyrinth  leak  rate. 
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Accordingly,  Robbins  determined  accurately  the  flow  rate 
into  the  stage  and  the  leak  rate  through  the  labyrinth  seal 
(See  Pig.  1)  for  all  operating  conditions.  The  results 
then  obtained  for  turbine  C were  reported  fully  in  Ref.  2. 
The  continued  presence  of  scatter  in  the  measured  loss 
coefficients  was  reported  and  a smoothing  technique  to 
eliminate  the  scatter  was  suggested. 

Before  the  Turbine  Test  Rig  could  be  used  to  measure 
the  effects  of  varying  parameters » 

1.  The  scatter  in  the  loss  coefficients  had  to  be 
eliminated. 

2.  The  overall  accuracy  of  the  performance  results 
evaluated  from  the  rig  measurements  had  to  be 
verified  in  some  way. 

The  resolution  of  these  problems  was  the  goal  of  the 
present  work  and  is  the  subject  of  this  report.  First, 
a satisfactory  method  was  found  for  smoothing  the  loss 
coefficients.  The  method  is  described  in  Section  III. 

The  approach  taken  to  verify  the  performance  of  the 
rig  was  to  first  devise  an  analysis  which  predicted  the 
performance  of  the  turbine  in  terms  of  unknown  loss  coeffic- 
ients, and  then  to  show  that  the  measured  loss  coefficients 
were  consistent  with  the  predicted  behaviour. 

A description  of  the  Turbine  Test  Rig  is  given  in 
Section  II.  The  analysis  of  the  behaviour  of  the  turbine, 
involving  a comparison  of  an  analytical  prediction  with 
the  results  of  a short  test  program,  is  described  in 
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Section  IV,  Details  of  the  analysis  and  the  computer 
program  are  given  in  Appendix  a. 
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II.  TURBINE  TEST  RIG  INSTALLATION 


A.  DESCRIPTION 

The  test  installation  consists  of  three  major  compon- 
ents! an  Allis-Chalmers  twelve  stage  axial  flow  compressor , 
an  exhauster  assembly,  and  the  turbine  test  rig  (TTR) 
itself. 

The  compressor  is  the  source  of  driving  air  for  the 
TTR  and  for  the  exhauster  assembly.  Fig.  2 shows  the 
piping  arrangement.  Turbine  air  passes  through  the  first 
settling  tank  into  an  eight-inch  pipe  containing  a flow 
nozzle,  into  the  second  settling  tank  and  into  the  turbine. 

Fig.  3 shows  the  plenum,  the  floating  stator  assembly, 
the  rotor,  and  the  dynamometer  (Ref.  1).  Pressure  ratios 
of  6«1  can  be  achieved  when  the  system  is  hooded.  The 
hood  was  needed  to  achieve  high  pressure  ratios  in  the 
tests  reported  here.  Fig.  4 shows  the  turbine  blading 
of  the  stator  and  rotor.  Ref.  3 contains  detailed  descrip- 
tions of  the  test  rig  hardware. 

The  floating  stator  assembly  shown  in  Fig.  3 permits 
measurements  of  the  axial  force  and  the  torque  on  the 
assembly.  Axial  and  rotational  movements  are  constrained 
by  calibrated  force  transducers  that  are  heat  insensitive. 
These  measurements,  together  with  wall  static  pressure 
measurements,  allow  the  determination  of  the  average 
axial  and  tangential  velocity  components  at  the  stator  exit. 

In  this  report  one  configuration  designated  Turbine  C 


tfas  tested,  the  geometry  of  which  is  shown  in  Fig.  5* 

Table  I describes  the  geometry  quantitatively.  The  stator 
blade  profile  is  shown  in  Fig.  6.  The  blades  of  the  stator 
generate  a converging-diverging  nozzle  shape.  Pressure 
measurements  were  taken  at  the  locations  shown  in  Fig.  6. 
The  pressures  necessary  to  the  analysis  of  the  stator  axial 
force  were  taken  at  the  locations  shown  in  Fig.  7 • 

B.  TEST  MEASUREMENTS  AND  ACCURACY 

1.  Mass  Flow  Rates 

Appendix  A of  Ref.  2 gives  a detailed  description 
of  the  method  used  to  determine  both  the  turbine  flow  rate 
and  the  labyrinth  seal  leak  rate. 

2.  Forces.  Torques,  Temperatures,  and  Pressure 
Ref.  3 and  Ref.  5 give  calibration  procedures 

for  the  TTR.  Identical  procedures  were  employed  here. 

Table  II  of  Ref.  1 gives  the  expected  accuracies  of  the 
measurements . 

C.  TESTING  AND  DATA  REDUCTION 

The  TTR  data  collection  system  is  described  in  Ref.  4. 
Appendix  D of  Ref.  1 gives  a detailed  explanation  of  the 
turbine  test  procedures.  Those  procedures  were  followed 
here  with  the  exception  that  a constant  RPM  was  held  and 
the  pressure  ratio  varied  over  the  desired  range.  The 
data  reduction  method  developed  in  Ref.  1 and  Ref.  2 was 
revised  as  described  in  Appendix  B. 
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III.  DATA  SMOOTHING  TECHNIQUE 

The  sensitivity  of  the  loss  coefficients  to  variations 
in  measured  quantities  was  shown  in  Ref.  2.  It  was  also 
stated  that  some  variation  is  unavoidable  since  measure- 
ments are  taken  over  a period  of  more  than  a minute.  It 
was  also  pointed  out  in  Ref . 2 that  the  parameter  most 
important  to  the  calculation  of  the  loss  coefficients  was 
Pj  (the  average  pressure  at  the  stator  exit).  is  an 
average  pressure  that  can  not  be  measured  directly.  It 
is  derived  as  described  in  Ref.  1 from  many  other  measure- 
ments. Significant  scatter  was  observed  in  the  variation 
of  P^  as  speed  was  varied  at  fixed  pressure  ratio.  However, 
it  was  found  that  the  hub  and  tip  pressures  (Ph  and  P^) 
measured  just  downstreajn  of  the  stator  varied  smoothly  over 
the  same  range.  Since  these  two  pressure  were  measured 
directly,  and  since  they  varied  smoothly,  it  was  assumed 
that  the  pressure  behind  the  stator  must  vary  smoothly  also. 
Consequently,  it  was  determined  that  a polynomial  curve  fit 
could  be  used  to  describe  the  variation  of  P^^  as  speed  was 
varied.  The  variation  of  Pj^  was  represented  as  a function 
of  Ph,  P^.,  and  the  isentropic  head  coefficient  (K^g) 
in  the  formi 

<r=^=K=fl  (1) 

Where  is  the  polynomial  coefficient. 

Then,  P^/P.^  (where  P^0  is  total  pressure  upstream  of 
the  stator)  was  computed  using  the  expression 
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Pig.  8 is  an  illustration  of  this  procedure,  which 
was  added  to  the  "bulk  process"  data  reduction  program. 

The  data  points  in  Fig.  8 were  taken  from  Runs  6 and  7 
in  Ref.  2.  The  two  Runs  were  for  the  same  conditions  but 
were  made  at  different  times.  It  can  be  seen  that  the 
trends  are  the  same  for  both  runs  but  the  scatter  is 
considerable.  The  scatter  is  the  cause  of  the  scatter  in  the 
calculated  loss  coefficients.  The  lines  on  Fig.  8 are  the 
polynomial  approximations  according  to  Eq.  (1)  for  the  two 
runs.  It  can  be  seen  that  the  polynomial  approximation 
averages  the  data  and  maintains  the  original  trend.  The 

S 

chosen  smoothing  function  was  a 1 degree  polynomial. 

Pig*  9 shows  the  stator  loss  coefficient  prior  to 
smoothing  for  the  points  in  Run  #7.  Fig*  10  shows  the 
same  data  after  smoothing. 

Similarly,  Fig.  11  shows  the  rotor  loss  coefficients 
for  Run  #7  prior  to  smoothing  and  Fig.  12  shows  the  same 
data  after  smoothing. 

It  can  be  seen  from  these  figures  that  the  scatter 
has  been  removed.  This  is  particularly  true  for  the  rotor 
loss  coefficients.  As  pointed  out  above,  it  is  believed 
that  the  observed  scatter  was  due  to  the  sensitivity  of 
the  loss  coefficients  to  small  changes  in  measured  quan- 
tities during  the  data  collection  process.  The  smoothing 
technique  removes  the  random  variations  recorded  during 
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data  collection  and  results  in  a much  more  realistic 
representation  of  the  variation  in  the  losses. 

The  smoothing  technique  was  incorporated  into  the 
data  reduction  program  as  described  in  Appendix  B. 
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IV.  ANALYSIS  OF  TURBINE  PERFORMANCE 

A.  APPROACH 

In  order  to  determine  if  the  performance  evaluated 

| • 

from  the  rig  measurements  was  accurate  an  analysis  to 
predict  the  behaviour  of  the  test  turbine  was  carried  out 
and  programmed  in  BASIC  language. 

A performance  test  was  then  conducted  in  a particular 
way  in  order  to  provide  a comparison  of  the  measured  with 
the  predicted  behaviour. 

B.  ANALYTICAL  PREDICTIONS 

The  flow  through  the  turbine  was  analysed  using  a 
pseudo-1  Dimensional  compressible  approach.  The  analysis 
is  described  in  detail  in  Appendix  A,  together  with  the 
computer  program  which  was  used  to  obtain  predictions  of 
the  turbine  performance.  One  of  the  inputs  which  the 
program  requires  is  the  rotor  passage  loss  coefficient  at 

: i 

zero  incidence.  Using  the  method  given  by  Vavra  in  Ref.  6 
it  was  determined  that  the  rotor  loss  coefficient  should 
have  the  value  .2514.  However,  results  of  previous  tests 
of  the  turbine  indicated  that  the  rotor  loss  coefficient 
was  rarely  as  high  as  .2514  and  could  be  as  low  as  .1. 
Therefore,  the  performance  of  the  turbine  was  analysed 
using  rotor  loss  coefficients  of  .2514  and  .1.  The  pre- 
diction program  was  run  for  both  15.000  and  18,000  RPM  with 
an  assumed  rotor  loss  coefficient  of  .2514,  and  for  18,000 
RPM  with  an  assumed  rotor  loss  coefficient  of  .1 
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The  above  parameters  were  chosen  to  obtain  a prediction 
of  the  turbine  performance  in  a range  in  which  experimental 
data  could  be  obtained.  In  particular,  the  analysis  could 
be  used  to  predict  the  pressure  ratios  at  which  choking 
would  occur  in  the  rotor  as  well  as  in  the  stator.  The 
pressure  ratio  at  which  choking  occurred  in  the  rotor  could 
then  be  established  experimentally  in  a test  conducted  at 
fixed  speed.  An  examination  of  the  choking  condition  was 
considered  to  be  a first  test  of  the  performance  analysis. 

The  results  of  the  analysis  for  the  parmameters  given 
above  are  shown  in  Pig.  13,  Pig.  14,  and  Pig.  15.  What  is 
shown  in  the  figures  is  the  map  of  the  range  of  values 
which  unknown  parameters  in  the  analysis  can  have,  that 
leads  to  a solution. 

Figures  16-19  show  predicted  performance  parameters 
for  the  case  of  18,000  RPM  and  assumed  rotor  loss  of  .1. 
These  results  will  be  discussed  in  conjunction  with  the 
results  of  the  turbine  test  run. 

C.  EXPERIMENTAL  TURBINE  TEST 

In  order  to  examine  the  occurrence  of  rotor  choking 
the  turbine  was  run  at  constant  RPM  and  the  pressure  ratio 
across  the  stage  was  increased  in  increments  by  lowering 
the  back  pressure.  The  point  at  which  the  horsepower  ceased 
to  increase  for  an  increase  in  pressure  ratio  was  examined 
to  determine  the  choking  point.  At  the  condition  where 
the  flow  reaches  a Mach  number  of  unity  at  the  exit  of  the 
rotor,  the  power  produced  by  the  turbine  can  not  be  changed 
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by  altering  the  downstream  pressure. 

Tests  were  conducted  in  this  manner  at  15,000  and 
18,000  RPM.  The  controlled  parameters  for  the  tests  are 
given  in  Table  II.  The  reduced  data  is  given  in  Table  III. 
The  referred  horsepower  is  shown  plotted  versus  the  pressure 
ratio  for  the  15,000  RPM  run  in  Fig.  20  and  for  the  18,000 
RPM  run  in  Fig.  21. 

The  results  are  discussed  in  the  next  section. 
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D. 


COMPARISON  OF  ACTUAL  AND  PREDICTED  PERFORMANCE 
Figures  13 » 14,  and  15  are  maps  of  possible  solutions 


for  the  flow  through  the  turbine,  when  the  parameters  which 
are  unknown  sure  allowed  to  vary.  It  is  noted  first  that 
the  predicted  range  of  solutions  extends  to  pressure 
ratios  telow  the  predicted  choking  line.  The 

explanation  for  this  apparent  anomaly  is  that  the  program 
calculates  all  possible  solutions,  and  for  any  point  below 
the  line  a particular  combination  of  losses  and  blockage 
factors  existed  which  would  allow  a solution  at  that  point 
without  choking.  It  should  be  pointed  out  that  any  point  on 
the  plot  cam  be  brought  to  the  choking  line  by  reducing 
the  stator  loss  by  a very  small  amount.  This  is  the 
procedure  that  was  followed  to  get  the  range  of  values 
at  choking  shown  in  Figures  16-19. 

Fig.  13  is  the  map  produced  by  the  program  for  RPM 
equal  15*000  and  for  an  assumed  rotor  passage  loss  coeffic- 
ient equal  to  .2514.  As  can  be  seen  from  the  figure  the 
pressure  ratio,  (Pg/P.^),  at  choking  was  about  .28,  cor- 
responding to  a stage  pressure  ratio  (P^/Pg)  °~  3*57. 

Fig.  14  is  the  map  for  RPM=18,000.  Note  that  the  predicted 
stage  pressure  ratio  for  choking  is  again  about  3.57* 

Fig.  20  and  Fig.  21  show  the  variation  of  the  referred 
horsepower  vs.  the  pressure  ratio  which  were  measured  in 
turbine  tests  at  15,000  and  18,000  RPM  respectively.  It 
can  be  seen  that  the  referred  horsepower  did  not  become 
independent  of  the  back  pressure  at  the  pressure  ratios 
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achieved  in  these  tests.  It  was  concluded  that  the  stage 
did  not  choke  at  the  predicted  pressure  ratio  of  3*57. 

It  was  noted  however  that  at  the  highest  pressure  ratios 
obtained  at  18,000  RPM,  the  slope  of  the  horsepower  curve 
was  becoming  smaller. 

The  velocity  diagrams  in  Figures  23-31  are  also  consis- 
tent with  the  argument  that  the  rotor  did  not  choke  during 
the  turbine  tests.  The  velocity  from  the  rotor,  V2, 
behaves  smoothly  and  in  a predictable  manner  throughout 
the  pressure  range  at  both  test  speeds.  This  might  not 
be  expected  if,  following  choking,  shock  waves  appeared 
downstream  of  the  rotor  exit  plane. 

Note  the  values  of  rotor  loss  coefficient  given  for 
the  test  results  in  Table  III.  With  the  exception  of 
point  #10,  which  was  considered  to  be  in  error,  they  were 
all  smaller  than  the  value  (.251*0  assumed  in  the  first 
performance  calculations.  As  explained  in  Appendix  A, 

' the  value  of  the  rotor  passage  loss  coefficient  used  in 

the  program  is  the  smallest  value  that  can  be  calculated 
for  the  overall  rotor  loss  coefficient.  Therefore,  if 
the  experimental  results  were  accurate,  the  computer  program 
could  not  predict  the  correct  performance  since  it  could 
never  calculate  a rotor  loss  less  than  the  input  value  of 
the  passage  loss,  which  was  .251**-.  It  was  as  a consequence 
of  this  observation  that  the  performance  was  re-calculated 
using  a rotor  passage  loss  coefficient  equal  to  .1. 

Figures  16-19  show  the  choking  behaviour  of  the 
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turbine  predicted  using  the  computer  program  for  an  assumed 
rotor  passage  loss  of  .1.  In  these  figures,  the  data  for 
test  point  #6  is  also  shown.  In  this  case  it  can  be  seen 
that  the  predicted  pressure  ratio  (Pg/^to^  at  cholcin6  is 
about  .245,  corresponding  to  a stage  pressure  ratio  (P^0/P2) 
of  about  4.08.  The  latter  is  slightly  higher  than  the 
maximum  pressure  ratio  that  was  attained  in  the  turbine 
tests.  (At  18,000  RPM  the  highest  attainable  pressure 
ratio  was  3*97.  due  to  the  characteristics  of  the  dynamom- 
eter. ) 

Pig.  17  is  the  range  of  efficiencies  predicted  for  the 
choked  condition.  It  is  noted  that  the  efficiency  measured 
at  test  point  #6  was  reasonably  close  to  the  predicted 
maximum  efficiency. 

Pig.  18  shows  the  range  of  stator  losses  for  which 
solutions  existed,  in  comparison  with  the  value  of  stator 
loss  measured  at  test  point  #6.  It  can  be  seen  that  the 
measured  value  intersects  the  predicted  range  of  possible 
solutions.  Pig.  19  shows  the  predicted  range  of  the  rotor 
coefficient.  Again,  it  can  be  seen  that  the  measured  value 
at  test  point  #6  overlaps  the  predicted  range  of  possible 
solutions. 

Pig.  16  shows  the  predicted  range  of  horsepower 
compared  with  the  horsepower  measured  at  test  point  #6. 

The  measured  horsepower  was  slightly  greater  than  the  max- 
imum value  which  was  predicted. 
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E.  DISCUSSION 

The  results  shown  in  Figures  16-19  illustrate  the 
uncertainty  in  the  prediction  of  the  performance  of  the 
turbine  when  only  the  rotor  passage  loss  and  blockage 
factor  are  known.  It  is  recalled  that  the  analysis 
satisfies  only  continuity  through  the  stage,  and  values  of 
an  additional  loss  coefficient  and  an  additional  blockage 
factor  must  be  established  to  obtain  a unique  solution. 

It  is  the  aerodynamic  shaping  of  the  surfaces  which  deter- 
mines these  factors.  Figures  16-19  show  the  possible 
range  of  performance  that  results  simply  from  the  areas 
of  the  passages  and  the  blade  angles. 

The  irregular  shapes  of  the  bounds  on  the  possible 
solutions  in  Figures  16-19  are  interesting.  There  is  no 
obvious  explanation  for  the  reduced  range  of  solutions  near 

kbl  = *81, 

The  results  obtained  in  the  present  work  suggest 
that  the  performance  of  the  turbine  as  measured  in  the 
turbine  test  rig  is  reliable.  It  had  been  thought  prev- 
iously that  the  rotor  loss  measurements  were  too  low. 

Here,  an  analysis  was  carried  out  to  predict  the  performance 
of  the  turbine,  and  the  predicted  results  have  shown  very 
good  agreement  with  the  performance  measurements  in  the  test 
rig. 

Whether  or  not  the  rotor  was  choking  had  also  been  a 
question  in  the  past.  In  the  present  work  both  the  analysis 
and  experimental  tests  have  shown  that  the  rotor  was  not 
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choked  in  any  test  to  date.  The  computer  program  with  an 
assumed  rotor  loss  of  .1  has  shown  that  the  test  rig  should 
not  be  choked  at  any  pressure  ratio  below  about  4.08.  Test 
data  has  shown  that  the  rotor  was  not  choked  at  a pressure 
ratio  just  slightly  below  4.0. 

It  is  also  of  interest  to  note  that  for  a rotor 
passage  loss  coefficient  of  .2514  the  predicted  pressure 
ratio  was  shown  to  be  independent  of  RPM,  At  both  15.000 
and  18,000  RPM  the  choking  pressure  ratio  (Pto/P2)  WaS 
calculated  to  be  about  3.57.  More  results  are  needed  to 
confirm  that  the  choking  pressure  ratio  is  indeed  indepen- 
dent of  speed. 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


1.  The  rotor  of  Turbine  C was  not  choked  at  any  pressure 
ratio  tested  so  far.  This  conclusion  is  based  on  the 
prediction  of  the  computer  program  and  also  on  the 
turbine  test  results. 

2.  The  loss  coefficients  measured  in  the  TTR  and  smoothed 
as  described  in  this  report  can  be  aocepted  as  being 
truly  representative  of  the  losses  in  the  separate 
blade  rows  of  the  turbine. 

3.  In  particular,  the  results  from  the  computer  program 
suggest  that  the  magnitude  of  the  measured  rotor  loss 
coefficients  is  probably  correct.  The  method  used 

by  Vavra  in  Ref.  6 predicts  much  larger  values  of  loss 
coefficients  for  the  present  rotor  geometry  than  those 
which  were  measured.  When  the  measured  value  of  the 
rotor  loss  coefficient  was  entered  into  the  computer 
• program,  there  was  good  agreement  between  all  the 
calculated  and  the  measured  turbine  performance 
parameters.  When  the  higher  (calculated)  value  of  the 
rotor  loss  was  entered,  there  was  a pronounced  dis- 
agreement between  the  calculated  and  the  measured 
performance . 

4.  The  computer  program  should  be  used  (and  progressively 
developed)  in  conjunction  with  all  tests  carried  out 
in  the  turbine  test  rig.  Most  importantly,  a test 
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should  be  conducted  to  determine  the  rotor  choking 
condition  experimentally.  This  will  require  the  purchase 
of  a water-brake  dynamometer  to  extend  the  power-speed 
range  of  the  test  rig. 

5.  Since  the  accuracy  of  the  test  rig  measurements  is 

no  longer  in  doubt,  experiments  to  determine  the  effect 
of  parameter  changes  (e.g.  axial  and  tip  clearances) 
on  turbine  performance  can  go  ahead. 


TABLE  I 


r 


TURBINE  GEQMETRT 


TURBINE 

STATOR 

ROTOR 

C 

1 

1 

DESCRIPTION 

NUMBER 

OP 

BLADES 

""(DO 

h(ia.) 

STATOR  1 

CONVERGING 

DIVERGING 

NOZZLES 

2.9058 

31 

4.184 

0.5775 

ROTOR  1 

CIRCULAR 

ARC 

7.119 

60 

R^-4.250 

h.,-0.732 
h2  *0.8475 

TABLE  II 

PARAMETERS  FOR  TURBINE  TEST 


POINT  # 
1 

2 

3 

4 

5 

6 


RPM 


15,000 

2.01 

15,000 

2.42 

15,000 

2.93 

15.000 

3.49 

15,000 

3.68 

18,000 

3-94 

18,000 

3-47 

18,000 

2.98 

18,000 

2.52 

18,000 

2.02 

TABLE 

Ill  TURBINE  TEST 

POINT  #1 

PRESSURE  RATIO 

= 2.011859338 

REF  FLOW  RATE 

= 1.020000000 

REF  R P Pi 

= 14262.15617 

REF  ROTOR 

MOMENT 

= 162.5192926 

REF  HP 

= 23.19933152 

EFF.  T-S 

- 0.712728663 

STATOR  LO: 

33  THEOR. 

= 0.107106795 

STATOR  LO- 

3S  COEFF. 

- 0.220295023 

R 0 TOR  L 0 S : 

3 COEFF. 

= 0.198187353 

ROTOR  LOS: 

3 THEOR. 

= 0.386101836 

I SEN.  HEAD  COEFF. 

= 4.036270755 

TH.  DEG.  OF  REACT I 

ON  = 6.025306644 

ACT.  DEG. 

OF  REACT I ON=-0 . 084678 1 82 

• VELOCITY  TRIANGLE 

DATA 

Vi  = 973 

“?  •“  r* 

1 O d !' 

V2  = 264.23594 

Vfi 1 = 226 

65467 

VA2=  155.42281 

Vi  j 1 — 946 

93688 

VU2=  132.49932 

ALPHA  1= 

76. 539  8 

7 ALPHA  2=  4 

BETA  1 = 

60. 4441 

8 BETA  2 = -6 

POTNT  #2 

P R b S S U RE  RATI U 

= 2.42052738 3 

REF  FLOW  RATE 

= 1 . 820000000 

REF  RPM 

- 14309.0588 

REF  ROTOR 

MOMENT 

= 132.334797 

REF  HP 

= 30.04433385 

EFF.  T-S 

= 0.748970958 

S TAT  0 R L 0 : 

;S  THEOR. 

= 0.107100795 

S T A T 0 R L 0 

;S  COEFF. 

= 6. 183632754 

ROTOR  LOS: 

: COEFF. 

= y • 1 y •-*  53 1*  j!  5 4 

ROTOR  LOS: 

i 

: THEOR. 

= 0.455681440 

ISEN.  HEAD  COEFF. 

= 4 .'94 17653 12 

TH.  DEG.  OF  REACTION  =-3. 76562E-63 

-H 


ACT.  DEG.  OF  REACT  I ON=-0. 272796901 

VELOCITY  TRIANGLE  DATA 

VI  - 1126.06784  V2  = 184.36934 

VA1=  265.46949  VA2=  178.98521 

VU1=  1094.32844  VIJ2=  44.23666 


ALPHA  1= 
BETA  1 = 


76 . 36422 
64.62313 


ALPHA 
BETA  c 


13.88078 
•70. 79960 
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TABLE  III  TURBINE  TEST  RESULTS  (CONT. ) 

POINT  #3 


p p t -I 

SURE 

RATIO 

93 

•~i  c,'  -?  i*-,  f a 

<1  j i o f 4 

REF 

FLOW 

RATE 

= 

1. 

02 

0000000 

REF 

RF'M 

= 

14 

29 

7.06424 

REF 

ROTOR 

MOMENT 

If 

i . 

7256537 

REF 

HP 

- 

3f 

. 6 

368547 

EFF. 

T-S 

= 

0. 

76 

9522364 

5TAT 

OR  LO 

SS  THEOR. 

0. 

10 

710O795 

STAT 

OR  LO 

SS  COEFF. 

= 

0. 

17 

2065787 

ROTOR  L 0 3 

S COEFF. 

st 

0 . 

09 

3520223 

ROTO 

P LOS 

S THEOR. 

2 

0. 

43 

179246 

I SEN 

. HEAD  COEFF. 

5. 

3 3 2 9 4 2 9 

TH. 

DEG. 

OF  REACT  I C 

N 

=- 

0 . 

0312391 

ACT. 

DEG. 

OF  REACT  I 

ON 

=- 

0. 

2417955 

VELOCITY  TRIANGLE  DATA 


VI  = 1243.51552 

I|.' 

= 1 3 ■ 3 • J 7*  t*  • J 

VA1=  314.71395 

VA2 

= 206.33375 

VU1=  1203. 19324 

V LI  2 

= -71.64356 

ALPHA  1=  75.3996 

9 

ALPHA  2=  -19. 18090 

B E T A 1 = 64.49  9 3 

i 

BETA  2 = -71.73607 

POINT  #4 

PRESSURE  RATIO 

-= 

3. 49095 1547 

FEF  FLOW  RATE 

- 

1 . 020008000 

REF  RF'M 

14290.27554 

REF  ROTOR  MOMENT 

SS 

131 . 5713712 

REF  HP 

2 

41. 16922137 

EFF.  T-S 

2 

O. 753539626 

| 

STATOR  LOSS  THEOR. 

2 

8. 107108795 

STATOR  LOSS  COEFF. 

= 

0. 158063259 

ROTOR  LOSS  COEFF. 

2 

O. 1 34491135 

ROTOR  LOSS  THEOR. 

2 

0 . 427937152 

I 

j 

I SEN.  HEAD  COEFF. 

- 6.783669585 

TH.  DEG.  OF  REACTION 

=-0.014684530 

ACT.  DEG.  OF  REACT I0N--0. £13834738 


VELOCITY  TRIANGLE  DATA 


TABLE  III  TURBINE  TEST  RESULTS  (CONT.) 

POINT  #5 


PRESSURE  RATIO 

= 3.677359655 

REF  FLOW  RATE 

- 1 . 02000O000 

REF  ft PM 

— 14292.77432 

REF  ROTOR  MOMENT 

= 186.8682436 

REF  HP 

= 42.3776356 

EFF.  T-S 

s Q • r 5 ci 8 0o64  »• 

STATOR  LOSS  THEOR 

. = 0.187100795 

STATOR  LOSS  CO EFF 

. - 0.140979499 

ROTOR  LOSS  COEFF. 

= 0.192867968 

* 

ROTOR  LOSS  THEOR. 

= 0.435423887 

j 

I SEN.  HEAD  COEFF. 

= 6.958584953 

TH.  DEG.  OF  REACTION  = 3.90802E-03 

ACT.  DEG.  OF  REAC 

T 1 0 N = — 0 . 2 2 3 6 9 6 7 5 9 

VELOCITY  TRIANGLE 

DATA 

VI  = 1362.18537 

V2  = 301,42998 

V A 1 = 361.18549 

VA2=  250.00963 

VU1=  1313.45034 

VU2=  -168.39007 

alpha  i = 74,627 

53  ALPHA  2 = -33.96162 

BfcTfi  1 = 64  ^ 69083  BETA  2 = -71.82123 


POINT  #6 


PRESSURE  RATIC 

3 . 938522427 

REF  FLOW  RATE 

1 . 028000000 

REF  RPM 

17091 . 82469 

1 

REF  ROTOR  MOMENT 

169, 2232867 

I 

REF  HP 

45 . 89 1 66403 

EFF.  T-S 

0.784438361 

STATOR  LOSS  THEOR.  = 

O.  187180795 

I 

STATOR  LOSS  COEFF.  = 

0. 141347589 

ROTOR  LOSS  COEFF.  - 

8. 187722748 

ROTOR  LOSS  THEOR.  = 

0.421639292 

I SEN.  HEAD  COEFF. 

= 5.050982877 

TH.  DEG.  OF  REACTION 

= 0.184266156 

ACT.  DEG.  OF  REACTION®  8. 175212337 


VELOCITY  TRIANGLE  DATA 

VI  = 1318.93669  V2  = 269.74423 

VA1 = 334.28819  VA2=  260.54759 

VU1=  1274.94805  VIJ2=  -69.83484 

ALPHA  1=  75.30782  ALPHA  2=  -15.80438 

BETA  1 = 61.53822  BETA  2 = -70.56818 
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TABLE  III  TURBINE  TEST  RESULTS  (CONT. ) 


POINT  #7 

PRESSURE 
F'EF  FLOW 


RATIO 

RATE 


3.467558140 
1 . 02000O000 


REF  RP 

M 

= 

1 

7039.50315 

REF  ROTOR 

MOMENT 

= 

1 

55. 9660642 

REF  HP 

42. 16691707 

EFF.  T 

— s 

0 

. 734673443 

STATOR 

LOS 

S THEOR 

■ “* 

0 

. 107100795 

STATOR 

LOS 

S COEFF 

• ™ 

0 

. 165054929 

ROTOR 

LOSS 

COEFF. 

0 

. 128352524 

ROTOR 

LOSS 

THEOR. 

= 

0 

• 4 1 y y t»  y y y 5 

ISEN. 

HEAD 

COEFF. 

= 

4. 668481576 

T H ■ D E la  a U 

F REACT 

ION 

= 

0.O51715826 

ACT.  DEG. 

OF  REAL 

T ION= 

0.080815056 

VELOC I 

TV  TRIANGLE 

DATA 

VI  = 1 

280. 

67141 

.. 

234.80211 

V A 1 = 

326. 

93646 

VA 

■L.  * 

234.79915 

VU1=  1 

238. 

22418 

VU 

2 = 

-1 . 17934 

ALPHA 

1 = 

75.207 

20 

ALPHA  2= 

BETA  i 

= 

tl U » 69il 

43 

BETA  2 - ~ 7 1 

POINT  #8 


PRESSURE  RATIO 
REF  FLOW  RATE 
REF  RPM 

REF  ROTOR  MOMENT 
REF  HP 
EFF.  T-S 

S T A T U F'  L 0 S S THE  0 R . 
STATOR  LOSS  COEFF. 
ROTOR  LOSS  COEFF. 
ROTOR  LOSS  THEOR. 


— 2.974737 8 9 3 
= 1 . 020009080 
= 17058.7 0 9 3 3 
= 137.689044 

- 37 . 267 49796 
= 0 • i’’  > 4 7 r { 0 5 r 

= 0.107100795 
= 0.139140155 
= 0.033159985 
= 0.414274730 


I SEN.  HEAD  COEFF.  = 4.  169275243 
TH.  DEG.  OF  REACTION  = 0.038020949 
ACT.  DEG.  OF  REACTION-  0.096859972 


VELOC 

I TV 

TRIANGLE 

DATA 

VI  = 

1200 

. 60325 

V2  = 

219. 

93553 

V A 1 = 

298 

.71293 

VA2= 

203. 

69610 

y ! 1 1 s 

1 162 

, 84940 

VU2= 

69. 

4 0 y I--  -i' 

ALPHA 

1 = 

75.5931 

34  AL 

PHA 

2= 

BETA 

1 = 

59.51871  BE 

Tfi  2 

- 

13.39624 
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TABLE  III  TURBINE  TEST  RESULTS  (CONT. ) 


POINT  #9 


PRESSURE  RATIO 

— 

2.519408034 

REF  FLO W RATE 

- 

1 . 020000000 

REF  PPM 

- 

17170.861O9 

REF  ROTOR  MOMENT 

1 13. 78678 8 4 

REF  HP 

= 

31.00049275 

EFF.  T-S 

0.743362897 

STATOR  LOSS  THEOR. 

0. 1071O0795 

STATOR  L 0 S S C U E F F . 

■= 

0. 189630321 

R 0 TOR  LOSS  C 0 E F F . 

= 

O. 191292357 

ROTOR  LOSS  THEOR. 

= 

0. 413068844 

I SEN.  HEAD  COEFF. 

= 3.567656122 

TH.  DEG.  OF  REACTION 

= y . u 5 3 i b'  3 6 7 y 

ACT.  LEG.  OF  REACTION 

1=  0.  108307260 

VELOCITY  TRIANGLE 

DATA 

VI  = 1103.46913 

'•12 

= 249.75075 

VA1=  259.39382 

VA2 

= 183.84024 

VU 1 - 1072.54789 

V IJ  2 

— 1 3 9 ■ y j 0 3 y 

ALPHA  1=  76.4041 

a 

ALPHA  2=  42, 

BETA  1 - 57.9474 

Y 

BETA  2 = -69, 

POINT  #10 

PRESSURE  RATIO 

= 

2.013415708 

REF  FLOW  RATE 

= 

1 . 028008000 

REF  RF'M 

= 

17153.52871 

PEF  ROTOR  MOMENT 

82. 24688 8 4 5 

REF  HP 

= 

22.39168077 

EFF.  T-S 

= 

0. 688130968 

STATOR  LOSS  THEOR. 

ss 

8.  107180795 

S T A T 0 R L 0 S S COEFF. 

= 

8.1141 13589 

ROTOR  LOSS  COEFF. 

= 

0.489230888 

ROTOR  LOSS  THEOR. 

s 

0.351684738 

I SEN.  HEAL  COEFF. 

= 2 . 787522827 

TH.  DEG.  OF  REACTION 

= 8.114732223 

ACT.  DEG.  OF  REACTI0N=-1 . 15697E-03 

VELOCITY  TRIANGLE 

DATA 

VI  = 987.09080 

I I •“> 

345.21785 

VA1=  210.24041 

VA2=  153.85315 

VU 1 - 964 . 44140 

VU2*  309.03475 

ALPHA  1=  77.7023 

• i 

ALPHA  2=  63 

BETA  1 = 55.6133 

>c 

• J 

BETA  2 = -66 
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FIGURE  8 ILLUSTRATION  OF  SMOOTHING  USING 
FIRST  DEGREE  POLYNOMIAL 
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FIGURE  11  ROTOR  LOSS  VS.  ISENTROPIC  HEAD 
COEFFICIENT,  BEFORE  SMOOTHING 
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FIGURE  13  PREDICTED  PERFORMANCE  RANGE » 

RPM=15,000,  INPUT  ROTOR  LOSS=.25l4 


B.7B  B.72  B.7M  B.7E  B.7B  B.BB  B.B2  B.BH 


PIG.  14 

PERFORMANCE  MAP 


4? 


0.70  0.72  0.7H  0.7E  0.7B  0.B0 


FIGURE  15  PREDICTED  PERFORMANCE  RANGE, 

RPM=l8 , OOO,  INPUT  ROTOR  LOSS=.l 
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FIGURE  16  PREDICTED  REFERRED  HORSEPOWER  VS.  K 
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FIGURE  17  PREDICTED  EFFICIENCY  VS. 
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FIGURE  18  PREDICTED  STATOR  LOSS  VS.  K 
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FIGURE  19  PREDICTED  ROTOR  LOSS  VS.  Kfc, 


PREDICTED  ROTDR  LOSS  COEFFICIENT 
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FIGURE  20  REFERRED  HORSEPOWER  VS.  PRESSURE  RATIO 
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FIGURE  21  REFERRED  HORSEPOWER  VS.  ERESSURE  RATIO 
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APPENDIX  A 

TURBINE  PERFORMANCE  ANALYSIS 
A-l  INTRODUCTION 

This  appendix  gives  a detailed  description  of  the 
theory  and  analytical  technique  used  to  predict  the 
performance  of  a turbine  in  the  Turbine  Test  Rig.  Variables 
used  in  this  appendix  are  defined  as  they  are  introduced. 

Also  given  are  the  details  of  the  computer  program  used 
to  implement  the  analytical  method.  The  variables  used 

in  the  computer  program  are  listed  and  defined  in  Table  A-I. 

% 

A-2  ASSUMPTIONS 

The  following  assumptions  were  mades 

1.  The  stage  was  choked  in  the  stator. 

2.  In  solving  the  continuity  equation  between  the 
stator  and  rotor  the  solution  corresponding  to  the  lower 
static  pressure  was  taken. 

Assumption  (1)  was  made  on  the  basis  of  past  exper- 
ience with  the  test  rig.  Ref.  1 states  that  there  is 
supersonic  flow  in  the  stator  and  an  examination  of  the 
pressure  distribution  in  the  stator  passages  confirms  this 
statement.  Fig.  A-l  is  a plot  of  the  pressure  distri- 
bution through  the  stator  by  pressure  tap  number  (See  Fig.  6). 

It  can  be  seen  in  this  plot  that  there  was  a sharp  pressure 
rise  between  tap  #3  and  tap  #4,  which  indicates  the 
presence  of  a shock  in  the  divergent  section  of  the  passage. 

Also,  the  level  of  static  pressure  at  the  throat  taps 
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in  comparison  to  the  supply  pressure,  was  observed  to  be 
independent  of  downstream  pressure  for  all  operating  con- 
ditions. 

Assumption  (2)  limits  the  range  of  possible  solutions 
as  shown  in  Fig.  A-2.  This  figure  is  a plot  of  flow 
coefficient  (jZf)  vs.  pressure  ratio  (P^/P^Q)  for  various 
loss  coefficients  (z).  (The  flow  coefficient,  a non-dimen- 
sional mass  flux,  is  introduced  in  Section  A-2.)  There 
are  two  possible  pressures  for  a given  flow  coefficient 
and  given  loss  coefficient.  Solutions  were  limited  to 
the  left  side  of  the  line  connecting  the  points  of  maximum 

flow  function.  The  procedure  was  chosen  because  it  resulted 

f 

in  the  "supersonic  root"  for  the  flow  from  the  stator. 

Also,  conditions  corresponding  to  choking  were  sought 
for  the  rotor,  and  the  locus  of  locally  sonic  conditions 
is  as  shown  in  Fig.  A-3.  It  should  be  noted  that  because 
of  the  definitions  of  the  flow  coefficient  and  loss  coef- 
ficient, the  "choking"  condition  at  which  the  downstream 
pressure  has  no  effect  upstream  of  the  plane  in  question, 
does  not  correspond  to  the  maximum  mass  flux  from  given 
stagnation  conditions. 

The  analysis  is  pseudo-l-dimensional.  Blockage  factors 
are  introduced  with  the  physical  cross  sectional  area  to 
account  for  non-uniform  flow  conditions.  Loss  coefficients 
are  defined  on  the  basis  of  kinetic  energies.  The  analysis, 
which  requires  only  mass  flow  continuity  through  the  stage, 
is  divided  into  three  stages;  the  stator,  the  interblade 
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region  and  the  rotor. 

While  the  analysis  itself  is  general,  solutions  can 
only  be  obtained  by  assigning  values  to  chosen  parameters. 

In  the  solutions  reported  here,*  the  flow  angle  at  the  stator 
exit  (*j) , the  blockage  factors  for  the  stator  throat 
(k^g)  and  rotor  exit  plane  (kbR),  and  ro^or  passage 
loss  coefficient  (zR)  were  given  values.  The  value  of 
was  fixed  ato^  =75°  based  on  the  results  of  the  turbine 
tests  and  from  the  blading  geometry.  The  blockage  factors 
were  set  at  ^^=.965  and  ktR=.85.  The  rotor  passage  loss 
coefficient  (zR)  was  an  input  variable  which  was  set  at 
.251^  (theoretical  value)  or  .1  (lowest  value  obtained  in 
test  results).  The  "rotor  loss  coefficient"  output  from 
the  program,  deludes  all  losses  from  the  stator  exit 

plane  to  the  rotor  exit  plane.  z\~2  corresPonds  "to  the 
"rotor  loss  coefficient"  evaluated  from  test  rig  measure- 
ments! consequently,  z^_2^zR. 
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A-3  BASIC  RELATIONS 


The  sketch 
shows  a general 
adiabatic  process 
with  entropy  increase 
(losses)  on  a T-S 
diagram  on  which 
the  temperature 
scale  has  been  divided 

by  the  constant  stagnation  temperature  of  the  process.  A 
perfect  gas  is  assumed.  (0)  represents  the  initial  stag- 
nation state  and  ( ) represents  the  final  state. 

The  non-dimensional  velocity,  X,  is  defined  as 
X = -^2-,,  where  V^Q  =v^C^f^9is  the  "total"  or  "limiting" 
velocity,  and  the  loss  coefficient,  z,  is  defined  as  shown. 
Note  that  the  stagnation  velocity  is  constant  throughout 
the  process. 

The  "flow  coefficient",  is  defined  here  as 
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' U.  14,  A fck 

where  w is  the  flow  rate,  is  the  density  at  the  initial 
stagnation  temperature  and  pressure,  A is  the  flow  area  and 
is  the  blockage  factor.  The  flow  coefficient  defined 
in  this  way  is  a non-dimensional  flow  rate  per  unit  area,  or 
mass  flux,  referred  to  initial  stagnation  conditions. 

The  flow  rate  is  given  by 

fA  V 
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so  that  Eq.  A(l)  can  be  written 

which,  in  terms  of  the  non-dimensional  velocity  defined 
above,  becomes 


The  loss  coefficient  is  defined,  as  shown  in  the  above 
sketch,  as 

a = T-rTis  = / - -X*  Aj?3) 

rtt  -Tis  X% 

Rearranging  Eq.  A(3), 


/- 

7 1 • 


and  using  the  isentropic  relationship  between  pressure  and 
temperatures 


A(4) 


Using  Eq.  A(4),  Eq.  A(2)  becomes  a single  equation 
for  0 as  a function  of  pressure  ratio  P/P.^.  In  analysing 
the  turbine,  generally  the  flow  rate  is  known  so  that  the 
flow  coefficient  can  be  calculated  using  Eq.  A(l).  On 
specifying  a value  for  the  loss  coefficient, ' ( z) , the  pres- 
sure ratio  can  be  calculated  from  Eq>  A(2)  using  Eq.  A(4). 
An  iterative  technique  is  used  (Newton's  Method)  starting 
with  an  initial  estimate  of  the  pressure  ratio.  Fig.  A(2) 
shows  0 as  a function  of  P/P^  for  various  values  of  z, 
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i. 


in  graphical  form.  The  solution  to  the  left  of  the  maxima 
is  obtained  by  beginning  the  iteration  with  a small  value 

°f  VPto‘ 

The  analysis  of  the  turbine  begins  by  specifying  that 
the  flow  rate  is  set  by  the  choking  of  the  stator  nozzles. 
The  pressures  at  successive  stations  are  then  calculated 
in  turn  from  the  flow  coefficient  as  described  above. 

A-4  STATOR  EXIT  CONDITIONS 

The  flow  coefficient  at  the  throat  of  the  stator  is  the 
value  of  the  flow  coefficient  corresponding  to  locally 
sonic  conditions  at  the  minimum  area.  In  the  present  work, 
the  maximum  flow  coefficient  at  a loss  coefficient  equal  to 
.05  (defined  as  0*)  was  taken  at  the  stator  throat, 
throughout  the  calculations. 

Since 


0*  = 


te—r 

Jt»  Vt»  A$  kys 


A(5) 


where  Ag  and  kb3  are  the  area  and  blockage  factor  at  the 
throat  station,  the  flow  rate  being  constant  through  the 
machine  requires  that 


W = ft.  Vt.  As  = constant 


at  stations  ahead  of  the  rotor,  where  T^Q  is  constant;  this 

implies  that 

As  ^Ls  0*  = constant 
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The  value  of  0*  can  be  obtained  analytically  when  the 


loss  coefficient  is  known.  The  maximum  occurs  where 

AA.  = 0,  which  leads  to  the  solution  (denoting  values  at  the 

stator  throat  by  an  asterisk), 


(X*)2  = B -Jb2  - C 
where  B = » C = 

"*  p,/ptc  • 
where  D 


A(6) 


A(7) 


Then  0*  is  given  by  Eq.  A(2). 

With  0*  known,  continuity  of  flow  rate  requires  that 
at  the  stator  exit, 


0 ,0.  Ain — 

1 *1  kj,J  C° * *» 


A(8) 


where  subcript  1 represents  the  stator  exit  plane  (station  1) 
and  •<,  is  the  flow  angle  at  thestator  exit. 

Knowing  0^ , P1/PtQ  can  be  found  by  iteration  as 
described  in  section  A-3»  Solutions  on  the  left  side  of 
Fig.  A-2  were  selected  by  beginning  the  iteration  with 
Pj/fc^  close  to  zero. 

After  obtaining  thevalue  for  then  X^^  can  be 


71 


I 


calculated  using 


ffrfl 

Then,  Pj/P.^  can  be  found  from 
pi/ptt  * 


A(9) 


A(10) 


and 


Ptl/Pto  = <Pl^to>  <Ptl/Pl) 


A(  11 ) 


A-5  INTERBLADE  SPACE 
The  interblade 
space,  from  Station 
(1)  to  Station  (i) , 
is  shown  in  the  sketch. 

Conditions  at  (i) 
are  calculated  in  a manner 


similar  to  those  at  Station  (1). 


k First,  by  continuity, 


The  value  of  c<i  is  obtained  from  the  condition  that  ang- 
ular momentum  is  conserved  in  the  interblade  space.  Then 


Xui  = “r7"X1  sine*' 


A(13) 


where  R^mean  radius  at  Station  1,  and  R^=mean  radius  at 


Station  i. 


Then 


= sin"1  ili 


A(l4) 


and,  also 


A(15) 


where  N is  the  RPM,  and  the  tangential  velocity  (Vu^) 
has  been  non-dimensionalized  as 


ui  v. 


A(16) 


Knowing  0^ , a value  for  P^/P.^  can  be  found  if  a 
value  of  the  interblade  loss  coefficient  (z^)  is  assumed. 
Then 


VPto  = <VPtl>  <Ptl/pto> 


A(17) 
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and 


A(18) 


With  the  absolute  flow  properties  established  at 
the  rotor  entrance,  the  relative  flow  conditions  can  be 
calculated.  First,  from  the  geometry  of  the  velocity 
diagram,  the  relative  flow  angle  ) is  given  by 


8,  - tan'1 

[Xi  cot  J 


A(19) 


where 


- 


A(  20) 


is  the  noh-dimensional  rotor  speed  at  radius  R^.  Then  the 
non-dimensional  relative  velocity  (Xwi)  is,  from  the  vel- 
ocity diagram,  given  by 


y _ cos  <Xl. 
wi  “ cos  Qi 


A(  21 ) 


and  the  temperature  , by 


A(  22) 
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The  equivalent  temperature  (Ref.  1),  is  given  by 


A(23) 


A(24) 


where 

Tgi  = equivalent  temperature  into  rotor  (See  Pig.  A-?) 
Ppi  = equivalent  pressure  into  the  rotor  (See  Fig.  A-?) 

A-6  ROTOR  CONDITIONS 

Using  the  values  calculated  thus  far  it  is  now  possible 
to  calculate  the  conditions  in  the  rotor.  Continuity 
requires  that 


where  a subscript  R denotes  the  exit  plane  in  the  rotor 
blading. 

It  should  now  be  noted  that  the  process  in  the  rotor 
frame  from  the  equivalent  conditions  at  (i)  to  the  exit 
of  the  rotor  is  entirely  similar  to  the  basic  process 
described  in  Section  A-3.  However,  the  equivalent  temper- 
ature takes  the  place  of  stagnation  temperature  and 
relative  velocity  replaces  velocity  in  the  given  equations. 


: 
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The  pressure  ratio  which  satisfies  the  flow  coefficient 
given  by  Eq.  A (25),  for  an  assumed  rotor  passage  loss 
coefficient  (zR)  is  Pg^Ei*  This  is  obtained  as  described 
in  Section  A-3.  The  corresponding  non-dimensional  velocity 
is  now 


A(  26) 


where 


t 

w2  • 


A (27) 


Hence 


Xw2  = Yw2, 


/lEl_ ' 
/**>  • 


A(  28 ) 


Prom  the  velocity  diagram, 


x|  = (r2  - Xw2  sin^)‘  + (Xw2  eos^) 


A(29) 


where  = 71°  (from  the  geometry  of  the  rotor),  and 
U2  = ^i*!^  • 

The  temperature  at  the  rotor  exit  is  given  by 


to 


(l-Zg) 


U„\ 

p2-)*  zr  b^A 
. Ey  J\to/ 


Tt2  ^2  2 

so  that  + X2 


A(30) 

A(31) 
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Next  the  loss  coefficient  that  includes  all  of  the 
losses  in  the  stage  from  Station  1 to  Station  2 can  be 
calculated.  This  coefficient  includes  all  of  the  losses  in 
the  inter-blade  region  as  well  as  those  in  the  rotor 
passage.  It  is  called  z^_2  here  in  fact  corresponds 

to  the  measured  rotor  loss  coefficient  in  the  TTR  since  all 
losses  aft  of  the  stator  in  the  TTR  are  included  in  the 
"rotor  loss  coefficient." 


. - , TE'T2 

*1-2  ' 1 Te-T2 


is 


A-7  STAGE  PERFORMANCE  PARAMETERS 


A(32) 


In  the  notation  used  for  the  turbine  test  rig  in 
Ref.  1 and  Ref.  2,  the  following  equations  determine  the 
turbine  performance  parameters  from  quantities  calculated 
in  the  above  steps t 


A(33) 
A(3^) 
A(  35 ) 
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A-8  THE  COMPUTER  PROGRAM 
A-8.1  Overview 


There  are  nine  degrees  of  freedom  in  the  com- 
puter program i 

1.  , the  flow  angle  out  of  the  stator 

2.  the  stator  exit  blockage  factor 

3.  k^,  the  interstage  blockage  factor 
k^R»  the  rotor  passage  blockage  factor 

5»  zg,  the  stator  loss  coefficient 

6.  z^,  the  interstage  loss  coefficient 

7.  zR,  the  rotor  passage  loss  coefficient 

8.  N,  turbine  RPM 

9*  Tt0,  total  temperature  upstream  of  the  stator 

o<,  is  input  on  the  basis  of  turbine  test  results. 

At  high  pressure  ratios  the  value  of  was  measured  to 
be  75°-?6°.  Therefore,  the  value  of  o<,  was  set  and  held 
at  75° • The  values  of  RPM  and  TtQ  are  thosB  conditions 
for  which  a solution  is  being  sought. 

All  the  other  parameters,  the  rotor  passage  loss 
coefficient  (zR)  and  blockage  factor  (kbR)  were  given 
chosen  values.  This  choice  was  made  because  the  incidence 
losses  were  included  in  the  interblade  calculation.  It 
was  thought  that  the  parameters  describing  the  flow  inside 
the  rotor  passages  could  be  held  fixed  as  the  speed 
was  allowed  to  vary. 

Consequently,  there  are  four  degrees  of  freedom  with 
which  the  program  works.  The  goal  for  the  program  is  to 
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find  those  solutions  which  result  in  flow  coefficients  in 
the  rotor  that  are  less  than,  or  equal  to,  the  "choking" 
value.  This  is  done  in  the  program  by  iteration  on  the 
remaining  four  degrees  of  freedom.  Fig.  A-5  is  a block 
diagram  showing  schematically  the  iteration  process.  The 
order  of  iteration  is  (from  the  most  frequent  to  the  least 
frequent),  z^,  zs,  k^,  and  k^.  The  loss  coefficients  for 
both  the  stator  and  interblade  always  begin  at  zero  and  they 
are  incremented  as  the  program  seeks  solutions.  A graphical 
representation  of  a typical  iteration  is  shown  in  Fig.  A-4. 

It  can  be  seen  that  the  program  can  find  a solution  for 
the  condition  where  the  loss  coefficient  is  equal  to  zero 
(denoted  as  (1)  on  Fig.  A-4).  On  the  next  iteration  the 
program  will  set  the  loss  coefficient  equal  to  .1  and  it 
can  be  seen  that  a solution  still  exists  at  point  (2).  On 
the  next  iteration,  however,  it  can  be  seen  that  a solution 
does  not  exist  for  a loss  coefficient  equal  to  .2.  In  this 
case  the  program  will  go  back  and  add  a smaller  increment 
to  the  previous  loss  coefficient  (.1). 

Each  time  an  added  increment  results  in  a loss  coef- 
ficient greater  than  the  maximum  allowable  the  program 
subtracts  that  increment  and  adds  a smaller  one.  In  this 
way  the  program  eventually  finds  the  maximum  loss  coef- 
ficient which  will  yield  a solution.  At  this  point, the  given 
flow  coefficient  is  at  the  maximum  point  for  the  calculated 
loss  coefficient  line.  At  this  point  "choking",  as  defined 
here,  has  occurred. 
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It  should  be  pointed  out.  that  this  "choking"  con- 
dition (the  point  of  maximum  flow  coefficient  for  any 
given  loss  coefficient)  is  not  the  point  corresponding 
to  locally  sonic  conditions.  In  Fig.  A-3  the  locus  of 
points  of  maximum  flow  coefficient  and  the  locus  of  points 
of  local  sonic  flow  are  both  shown.  The  point  of  locally 
sonic  flow  always  falls  to  the  left  of  the  "choking"  point 
(as  defined  in  this  work),  and  the  flow  coefficient  for 
sonic  flow  is  always  less  than  the  flow  coefficient  at 
the  defined  choking  point.  The  reason  is  that  the  point 
of  maximum  mass  flux  (maximum  flow  coefficient)  is  for 
given  upstream  stagnation  conditions  (see  next  section), 
Because  of  the  way  the  flow  coefficient  is  defined,  the 
maximum  for  a given  loss  coefficient  will  not  correspond 
to  local  sonic  conditions  arxcept  where  the  loss  coefficient 
is  equal  to  zero. 

The  difference  between  the  "choking"  flow  coefficient 
corresponding  to  the  maximum  and  the  flow  coefficient  at 
the  local  sonic  condition  is  about  10%  at  a loss  coefficient 
of  .3  (the  highest  rotor  loss  coefficient  predicted  by 
the  program).  The  difference  is  about  4%  at  a loss  coef- 
ficient equal  to  .25  (the  highest  stator  loss  predicted), 
and  insignificant  (less  than  1%)  at  a loss  coefficient  of 
about  .15. 
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A-8 . 2 Description  of  the  Program 

The  following  discussion  refers  to  Fig.  A-6, 
which  is  a detailed  flow  chart  of  the  program,  to 
Table  A-I , which  is  a list  of  program  variables  and  their 
definitions,  and  to  the  program  listing  given  in  Section  A-10. 

The  first  step  in  the  program  is  to  input  the  variables. 
aH,  is  input  in  line  280  (it  is  also  converted  to  radians 
in  that  line),  k^  in  line  286  and  k^  is  input 

in  line  284.  All  other  variables  are  input  as  requested 
by  the  program  (lines  20-110).  Note  that  P2  and  P^.Q  are 
input  also.  They  are  not  involved  in  the  calculations  and 
were  originally  input  as  the  conditions  for  which  a 
particular  solution  was  sought. 

The  next  step  is  to  go  to  a subroutine  to  calculate 
the  value  of  0*  (line  170).  0*  is  defined  as  the  value 

of  the  flow  coefficient  at  the  throat  of  the  stator  which 
sets  the  values  of  the  flow  coefficients  at  all  downstream 
stations.  0*  is  generated  in  the  subroutine  by  the  method 
covered  in  Section  A-4  assuming  a loss  coefficient  of  .05 
to  the  throat  of  the  stator  (line  160). 

Next  the  program  calculates  the  conditions  in  the 
stator  exit  plane.  This  is  done  by  first  calculating  the 
value  for  0^  (line  320).  As  explained  earlier,  the  value 
for  0^  determines  the  upper  limit  of  the  stator  loss 
coefficient. 

Using  the  method  described  in  Section  A-3>  the  program 
calculates  the  value  for  . Note  that  if  the  stator 
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loss  coefficient  is  large  enough,  so  that  a solution  does 
not  exist,  this  will  be  detected  in  the  decision,  "Pj/P.to<0 
or>l"  (lines  350  and  370).  If  "yes"  then  the  stator  loss 
(zD)  is  too  large  and  the  following  steps  are  takeni 

1.  The  last  increment  added  to  the  stator  loss  is 
subtracted  (line  420).  (Note,  the  initial  stator  loss 
increment  is  .1.) 

2.  The  existing  stator  loss  increment  is  multiplied 
by  some  factor  less  than  1 (line  430). 

3.  The  new,  smaller,  increment  is  then  added  back 
to  the  stator  loss  coefficient  (line  1570). 

If  stator  loss  increment  is  sufficiently  small 
(line  1575)  then  the  upper  limit  of  the  loss  coefficient 
has  been  reached  as  discussed  earlier.  At  this  point  no 
further  solutions  are  attainable  and  k^  is  incremented 
(line  2290).  Note  that  the  upper  limit  of  stator  loss 
will  not  normally  be  attained  until  several  iterations 
have  been  made  through  the  entire  program. 

Next  the  interblade  calculations  are  made.  The  sol- 
ution for  P1/P^0  found  for  the  stator  exit  plane  is  used 
to  make  the  calculations  for  0^  (line  770) • Again,  0. 
determines  the  upper  limit  on  the  interblade  loss.  In 
the  interblade  region,  however,  there  is  another  criterion 
that  must  be  met.  The  value  for  Q generated  within  the 
program  must  be  equal  to  ($g  =69°  and  is  fixed  by  the 
geometry).  This  criterion  is  met  by  iterating  the  inter- 
blade loss  coefficient  as  follows  1 On  every  pass  through 
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this  section  the  interblade  loss  is  initially  set  to  zero 
and  the  interblade  loss  increment  is  initially  set  to  .1 
(lines  720,  740).  The  value  for  P^/P.^  is  calculated  using 
the  same  procedure  as  described  earlier  for  P^/P^.  Again, 
the  value  of  P./fc^  is  checked  on  each  iteration  to  ensure 
that  it  falls  between  0 and  1 (line  780).  If  it  does  not 
then  the  interblade  loss  is  reduced  in  a method  analagous 
to  that  covered  for  the  stator  loss. 

After  generating  a value  for  P^/P.^  the  program  cal- 
culates (line  1208).  When  the  value  for  is  generated 
it  is  compared  to  ^ and  three  results  are  possible. 

If  is  equal  to  f • 05  then  convergence  has  occurred  and 
the  program  will  continue.  If  is  less  than  gg  (line  1250) 
then  the  value  for  interblade  loss  coefficient  is  increased 
by  the  current  value  of  interblade  increment  and  the  program 
goes  back  to  the  start  of  the  interblade  calculations  and 
begins  anew.  If  ^ is  greater  than  (line  1240)  the 
value  of  interblade  loss  is  decreased  in  the  same  way 
as  was  described  for  the  stator  loss. 

Two  results  are  possible  in  the  interblade  regiom 
Either  the  program  finds  the  value  of  the  interblade  loss 
that  gives  flow  angle  convergence  (^;  = ) , or  the  inter- 

blade loss  increment  becomes  sufficiently  small  (10~^)  to 
trigger  termination. 

Note  that  the  interblade  loss  increment  is  allowed  to 
become  very  small  (10-^)  (line  850).  This  is  because  the 
above  mentioned  angle  convergence  requirement  is  very 
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sensitive.  Note,  also,  that  when  the  interblade  loss 
increment  falls  to  10“^  and  the  iteration  process  in  the 
interblade  region  is  terminated,  the  stator  loss  coefficient 
is  incremented  (line  850).  Stator  loss  is  incremented 
because  a higher  value  of  stator  loss  may  result  in  con- 
ditions that  allow  convergence  in  the  interblade  region. 

If  not,  the  stator  loss  calculations  eventually  reach  the 
termination  criterion  mentioned  earlier. 

Once  angle  convergence  has  been  attained  in  the  inter- 
blade region  then  the  value  of  the  rotor  flow  coefficient 
can  be  calculated  (line  1470).  If  the  program  gets  to  this 
point  for  any  given  set  of  conditions  then  a solution 
exists  and  the  only  determination  to  be  made  is  whether 
the  resulting  rotor  flow  coefficient  is  greater  than  the 
"choking"  value  (line  1530).  If  "yes"  then  it  is  dis- 
allowed, since  a flow  coefficient  greater  than  the  "choking" 
value  cannot  occur  physically.  If  the  calculated  coef- 
’ ficient  is  less  than,  or  equal  to  the  "choking"  value  then 

the  results  are  printed,  the  stator  loss  is  incremented,  and 
the  calculations  start  again  at  the  stator  exit  plane. 

The  choking  value  of  the  flow  coefficient  for  the 
rotor  (0*R)  is  determined  in  exactly  the  same  way  as  the 
choking  value  for  the  stator  (0*),  with  the  assumed  value 
of  rotor  loss  used  as  the  loss  coefficient.  (At  present 
the  value  for  ^ is  calculated  separately  and  inserted  in 
line  121.  It  is  suggested  that  a "GOSUB  2140"  be  inserted 
at  this  point.  Then  the  input  value  for  rotor  loss 
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coefficient  can  be  changed  easily  and  the  "choking"  value 
for  the  rotor  will  be  calculated  in  the  program.)  If 
the  resulting  flow  coefficient  is  greater  than  "choking", 
the  program  increments  stator  loss  coefficient  and  returns 
to  the  stator  plane  calculations  to  begin  anew. 

Therefore,  in  all  cases  the  program  will  seek  the 
highest  value  of  stator  loss  coefficient  for  which  solutions 
exist.  At  some  point,  the  rotor  flow  coefficient  normally 
attains  a value  less  than  the  "choking"  value  and  all 
solutions  for  higher  stator,  loss  coefficients  will  be 
acceptable  solutions.  Or,  if  the  "choking"  condition  is 
never  reached,  the  program  will  attain  the  highest  possible 
stator  loss  in  attempting  to  achieve  solutions  less  than 
"choking". 

When  the  highest  value  of  the  stator  loss  is  reached, 
there  are  two  possible  avenues  in  the  program.  If  the 
program  was  generating  solutions  less  than  "choking"  at 
the  time  of  termination  in  the  stator  calculations,  then 
kbi  will  be  incremented,  stator  loss  will  be  set  back  to 
zero,  and  the  process  begun  anew  (line  2280).  If  the 
solutions  being  generated  at  the  time  of  termination  were 
greater  than  "choking",  then  the  upper  value  of  kbi  has 
been  reached  and  any  increase  in  k^  will  yield  no  further 
solutions.  In  this  case  k.^  is  incremented,  k^  is  set 
to  the  pre -determined  lower  limit,  and  stator  loss  is 
reset  to  zero  (line  2330). 

The  program  continues  in  the  above  manner  until  the 
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pre -determined  upper  limit  of  kfel  is  reached  (line  2335) » 
at  which  point  the  program  is  terminated. 

Since  the  21-MX  computer  has  only  6 digits  of  accuracy, 
round-off  error  has  been  a problem.  Where  Newton's  Method 
is  used  to  solve  for  pressure  ratio  the  computer  lacks 
sufficient  accuracy  to  converge  to  the  desired  accuracy 
under  certain  conditions.  When  this  happens,  the  program 
will  stay  inside  the  iteration  and  never  converge.  This 
problem  has  been  solved  by  putting  a counter  inside  each 
iteration  loop  (lines  540  and  980).  More  than  60  iterations 
is  treated  as  a non-convergent  condition,  and  the  approp- 
riate loss  coefficient  is  reduced  in  the  manner  covered 
previously. 

A-9  OPERATING  PROCEDURE 

The  procedures  for  operating  the  Hewlett-Packard 
21-MX  computer  will  not  be  covered  in  this  paper  since 
the  applicable  manuals  are  available  at  the  laboratory. 

The  computer  must  be  on,  with  the  RTE-B  operating  system 
"READY-. 

1.  Load  the  paper  tape  program  labeled  "TTR  11". 

2.  Edit  into  the  program  the  minimum  values  of  kbl 
(line  284)  and  line  (286).  Edit  the  desired  value  of 
0<t  into  line  280.  If  it  is  desired  to  get  a print-out 

of  every  solution  then  either  remove  lines  1724  and 
1726  or  place  "REM"  in  front  of  them.  If  these  two  lines 
are  left  in  the  program  then  a solution  will  be  printed 
out  only  if  the  calculated  pressure  ratio  (P2^to^  ^ess 
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than,  or  greater  than,  the  previously  calculated  lowest  or 
highest  pressure  ratio  respectively,  for  the  value  of 
then  being  used. 

3.  Type  "RUN",  "RETURN"  on  the  keyboard  and  the 
program  will  begin  execution  by  asking  for  inputs  1 

"INPUT  STATOR  LOSS  COEFFICIENT"  The  normal  input  is 
zero.  However,  if  some  specific  case  is  required  then  any 
desired  loss  coefficient  can  be  input. 

"INPUT  PTO"  Input  the  upstream  total  pressure. 

"INPUT  P2"  Input  the  hood  pressure. 

"INPUT  RPM"  Input  the  RPM  desired. 

"INPUT  TTO"  Input  the  total  temperature  into  the 
stage . 

This  is  all  that  is  required  to  operate  the  program. 

Depending  on  the  range  of  k^  requested  the  program  will 
take  from  30  minutes  to  36  hours  to  run  on  the  21-MX. 

It  is  important  to  realize  that  the  above  procedure 
will  produce  solutions  that  give  values  of  the  rotor 
flow  coefficient  less  than  "choking".  If  the  turbine 

! 

performance  at  the  "choking"  cbndition  is  desired  then 
after  the  computer  run  is  complete  a further  step  must  be 
taken.  It  is  necessary  to  take  the  values  of  k^  and  k^ 
for  which  a solution  is  desired  and  force  the  choking 
condition.  This  is  done  as  follows  1 

1.  Scan  the  output  results  at  the  desired  k^  and 
k^  and  locate  the  point  that  has  a value  of  rotor  flow 
coefficient  closest  to  the  choking  value  and  a non-zero 
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stator  loss  coefficient. 

2.  Re-start  the  program,  this  time  putting  in  a value 
of  the  stator  loss  coefficient  slightly  less  than  the  value 
printed  on  the  output. 

As  stator  loss  is  manually  decreased,  the  rotor  flow 
coefficient  will  increase  towards  the  choking  value.  When 
the  rotor  flow  coefficient  is  within  the  limits  specified 
in  line  2370  the  words  "CONVERGENT  CONDITIONS"  will  appear 
on  the  teletype  output  preceding  the  printed  results. 

If  the  stator  loss  coefficient  is  decreased  too  much, 
the  rotor  flow  coefficient  will  increase  to  a value  above 
choking.  In  this  case  no  ouput  will  be  printed  on  the  tele- 
type. However,  the  flow  coefficient  will  appear  on  the 
video  display.  In  this  case  increase  the  stator  loss 
slightly. 

It  has  been  found  in  this  work  that  the  extreme  values 
of  referred  horsepower,  stator  loss,  rotor  loss,  and 
efficiency  occur  at  the  lowest  and  highest  values  of  k^ 
for  each  value  of  k^.  With  this  in  mind  it  is  possible  to 
find  the  range  of  predicted  values  at  choking  by  forcing 
to  choking  only  the  lowest  value  of  kbi  (with  non- zero 
stator  loss)  and  the  highest  value  of  kbi  for  each  k^. 
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A-10  PROGRAM  LISTING 
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AD-A035  876 
UNCLASSIFIED 


NAVAL  POSTGRADUATE  SCHOOL  MONTEREY  CALIF  F/G  21/5 

ANALYTICAL  AND  EXPERIMENTAL  DETERMINATION  OF  THE  CHARACTERISTIC— ETC (U) 
DEC  76  B C BOATRIGHT 
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4080  LET  D=P0/29.92 
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r~ 

i 

TABLE  A-I 

j 

DEFINITION  OF  VARIABLES 

AO  4A*)  Area  at  the  Stator  Throat 
A1  ( A Area  at  the  Stator  Exit  Plane 

A8  Decision  Variable 

A9  Decision  Variable 

Cl  Decision  Variable 


D4 

E 

Decision  Variable 

EO 

(^T-S^  Total- to-Static 

Efficiency 

FO 

(0*)  Flow  Coefficient 

at  Stator  Throat 

FI 

(0^)  Flow  Coefficient 

at  Stator  Exit  Plane 

F2 

(#R)  Flow  Coefficient 

Through  Rotor  Passage 

F3 

(0^)  Flow  Coefficient 

Through  Interblade  Area 

Fk 

Input  Value  of  0 Based  on  Assumed  Rotor  Loss 
R 

F8 

Decision  Variable 

F9 

Decision  Variable 

GO 

(Gamma  - 1)/Gamma 

HO 

Horsepower 

HI 

Referred  Horsepower 

H8 

(?)  Rho 

H9 

Reynolds  Number 

J 
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J 


in)  Pi 


K1  k£g  Blockage  Factor  at  Stator  Throat 

K4  kfei  Interblade  Blockage  Factor 

K5  kbR  Blockage  Factor  in  the  Rotor 

K6  k^^  Blockage  Factor  at  the  Stator  Exit  Plane 

Ml  Rotor  Torque 

M2  Referred  Rotor  Torque 

N RPM 

NO  Referred  RPM 

01  (&i)  Flow  Angle  at  Stator  Exit  Plane 

03  (tfj  Flow  Angle  at  Rotor  Entrance  Plane 

04  (Akfci^  Increment  for  kfei 

p ptl/pt0 

po  Pto 

Pl  V*to 

pz  P/Pto 

p3  VPtl 

P4  ptl/ptc 

p5  Pi/Ptl 

p«  Vpto 

p7  pEi/pto 

P8  Decision  Variable 

P9  P2 

R (jt«)  Total  Density  at  Stagnation  Conditions 

R7  Relaxation  Parameter 

R8  Relaxation  Parameter 

R9  Predicted  Pg/P^ 
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S7-S9  Used  in  GOSUB  7000  to  Find  Pg/P^ 

TO  Tt0 

11  Tt2/Tto 

T2  T2/Tt0 

T3  VTto 

TEi/Tto 

T5  VTto 

T6  TEi/Tto 

U1  Dimensionless  Velocity  (See  Appendix  A) 

U2  TTg  Dimensionless  Velocity  (See  Appendix  A) 

U3  ^ Dimensionless  Velocity  (See  Appendix  A) 

VO  VtQ  Total  Velocity 

WO  w Flow  Rate 

W1  w*  Referred  Flow  Rate 

W9  Decision  Variable 

XO  Xwl  Non-Dipiensional  Relative  Velocity 

XI  Non-Dimensional  Velocity 

X2  X2  Non-Dimensional  Velooity 

X3  Xj^  Non-Dimensional  Velocity 

X4  Xui  Non-Dimensional  Velocity 

X5  Xwi  Non-Dimensional  Relative  Velocity 

X6  Xw2  Non-Dimensional  Relative  Velocity 

X8  Decision  Variable 

Y2  Non-Dimensional  Velocity 

Y6  Increment  for  fc^ 

Y7  Increment  for  k^ 

Y8-Y9  Counters 


i 
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zo  zg  Stator  Loss  Coefficient 

Z3  z ^ Interblade  Loss  Coefficient 

27  Z1 -2  ?re dieted  Rotor  Loss  Coefficient 
Z8  Increment  for  Stator  Loss  Coefficient 
Z9  Increment  for  Interblade  Loss  Coefficient 


101 


FIGURE 


FIGURE  A-3 


IGURE  A-4- 
OLUTION  PROCESS 


FIGURE  *-?  THERMODYNAMIC  PROCESS  OF  FLUID 


IN  AN  AXIAL  TURBINE  STAGE 
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APPENDIX  B 

TURBINE  TEST  RIG  (TTR)  DATA  REDUCTION  AND  PROCESSING 
B-l  INTRODUCTION 

This  appendix  describes  only  those  changes  that  have 
been  made  to  the  data  reduction  procedure  given  in  Ref.  2. 
The  program  numbers  and  the  functions  of  those  programs 
given  in  Ref.  2 have  not  changed.  The  channel  and  port 
assignments  for  data  collection  and  storage  have  not 
changed.  Variables  also  have  not  changed,  although 
additional  variables  were  defined  in  the  course  of  mod- 
ifying the  data  reduction  process.  Those  variables  that 
now  exist  in  addition  to  those  given  in  Ref.  2 are  given 
in  Table  B-I. 

Modifications  made  during  the  course  of  this  work  weret 

1.  Raw  data  can  now  be  read  directly  from  mass  memory. 
Previously  it  was  necessary  to  read  the  paper  tape  for  a 
given  point  in  order  to  reduce  it. 

2.  Provisions  have  been  written  into  the  programs 
for  smoothing  the  reduced  data. 

3.  All  of  the  parameters  used  in  calculating  the 

theoretical  loss  coefficients  that  were  previously 
entered  from  charts  are  now  in  polynomial  form  in  the 
program.  

4.  Storage  of  raw  data  is  now  a separate  program. 
Previously  it  was  necessary  to  run  the  reduction  sequence 
to  store  data. 
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B-2  DESCRIPTION  OP  PROGRAMS 


The  data  reduction  is  divided  into  ten  separate 
programs  with  an  additional  program  to  store  the  raw  data. 
Fig.  B-l  shows  the  contents  of  each  program  and  also  shows 
the  sequence  in  which  the  programs  are  chained.  Following 
is  a description  of  the  reduction  sequence  referred  to 
Fig.  B-l  which  will  be  followed  by  a description  of  how  to 
run  the  reduction  program. 

TTR  in  Fig.  B-l  is  used  only  to  store  the  raw  data. 
TTR  does  not  chain  to  any  other  program. 

The  data  reduction  process  begins  in  TTR1.  In  order 
to  smooth  the  reduced  data  as  discussed  in  Section  III  it 
is  necessary  to  have  n values  for  (pi“p|,ub^^Ptip”Phub^  * 
Therefore,  the  first  step  in  the  data  reduction  process  is 
to  go  from  TTR1  to  the  point  in  TTR2  where  the  values  of 
Pl*  Ptip’  3111(1  Phub’  have  1)6611  calculated.  On  completion, 
the  program  returns  to  TTR1  and  repeats  the  process  n 

times  until  the  n values  of  (pi-Phub^Ptip“Phub^  have 
been  calculated.  The  program  then  proceeds  to  TTR9  where 
a polynomial  curve  fit  for  the  n values  of 

( P1  “Phub W Pt ip~Phub ) is  66116rated.  Then  the  program 
returns  the  polynomial  coefficients  to  TTR1B  where  the 
reduction  process  begins  at  the  first  point  in  the  run. 
However,  when  the  point  in  TTR 2 is  reached  where  is 

calculated,  it  is  calculated  using  the  polynomial  as 
described  in  Section  III. 

The  rest  of  the  reduction  process  in  Fig.  B-l  is  the 
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sane  as  described  in  Ref.  2 with  the  exception  that  no 
further  keyboard  inputs  are  needed  until  the  raw  and 


reduced  data  is  tabulated. 


B-3  RAW  DATA  STORAGE 

The  procedure  for  storing  raw  data  on  the  mass  memory 
is  outlined  in  Ref.  2(  Appendix  B.  Sec.  B-3*  lines  10t19* 

Note  that  line  10  should  read,  "GET  ’TTR'".  Although  it  is 
now  a separate  program  the  procedures  for  storing  data 
have  not  changed. 

B-4  PROCEDURE  FOR  DATA  REDUCTION  PROGRAM 

1.  Key  in  GET  "TTR1B" 

2.  Press  "RUN  EXECUTE". 

3.  The  following  check  list  will  be  printed  on  the 
HP  9830  printer 1 

"PRIOR  TO  RUNNING  THIS  SEQUENCE  ENSURE  FOLLOWING! 

"THAT  TTR 2 LINE  590  HAS  PROPER  FACTOR  IN  ITj 

1.02  FOR  HOODED 
1.01  FOR  UNHOODED 

"IF  IT  IS  NF CESS ARY  TO  INPUT  ALPHA  1 THEN  CHANGE  TTR2 
"LINE  960  TO  'INPUT  A3' 

"IF  BLOCKAGE  FACTOR  OTHER  THAN  1,  CHANGE  TTR2-1060 
"TO  'INPUT  X7',  AND  PUT  SEMICOLON  AFTER  TTR2-1050 
"ENSURE  THAT  TTR1 , TTR1B  HAVE  THE  PROPER  FILENAME  IN  441 
"IF  IT  ISN'T  DESIRED  TO  STORE  REDUCED  DATA  CHANGE 
TTR6-580  TO  'G1=0' 

"ENSURE  TTR6-610,  63 0 HAVE  PROPER  FILENAME  FOR  REDUCED 
DATA 


"ENSURE  TTR7-220  HAS  PROPER  FILENAME  FOR  RAW  DATA 

"ENSURE  TTR8-220  HAS  PROPER  FILENAME  FOR  REDUCED  DATA 

"IF  IT  IS  DESIRED  TO  RUN  A SINGLE  POINT,  WITHOUT 
SMOOTHING,  INSERT  TTR2-1127  'GOTO  1200'.  OTHERWISE 
OMIT  TTR2-1127 

The  above  checklist  will  prepare  all  ten  programs.  The 
last  item  is  a provision  for  elimination  of  smoothing  if 
it  is  desired  to  run  a single  point.  If  smoothing  is  not 
desired  then  each  point  must  be  run  separately  with 
TTR2-1127  inserted. 

4.  Press  CONTINUE  EXECUTE 

5.  The  display  will  read  LOWEST,  HIGHEST  RECORD  # 

THIS  RUN.  Input  the  lowest  record  number  and  the 

highest  record  number  from  which  it  is  desired  to  read 
raw  data.  Note,  the  record  number  on  which  the  reduced 
data  will  be  stored  will  be  the  same  number  as  the  raw 
data  record  from  which  the  raw  data  was  read.  Make  sure 
the  filename  designated  for  the  storage  of  the  reduced 
data  can  accept  data  on  those  record  numbers  without 
writing  over  previous  information. 

6.  Next  the  display  will  read  PRINT  OUT  RESULTS? 

YES=1,  N0=0.  If  a 1 is  selected  here  then  at  the  conclusion 
of  data  reduction,  for  all  the  points  in  the  run,  the 
program  will  begin  to  print  all  results.  If  0 is  selected 
then  all  reduced  data  will  be  stored  on  the  mass  memory. 


7.  It  will  take  about  15  minutes  for  the  program  to 
read  the  raw  data  and  make  the  calculations  necessary  for 


smoothing  the  data.  When  this  is  oompleted  the  program  will 
start  reducing  each  point  and  printing  the  results. 
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8.  At  the  completion  of  data  reduction  for  all 
points  the  option  to  tabulate  the  data  is  available. 

ENTER  RECORD  #'St  LOWEST,  HIGHEST  will  appear  on  the  display. 
If  it  is  desired  to  tabulate  the  raw  data  then  enter  the  same 
record  numbers  that  were  entered  to  initiate  the  program. 

The  raw  data  will  then  be  tabulated. 

9.  Next  the  display  will  read  ENTER  LOWEST,  HIGHEST 
RECORD  NUMBER.  If  it  is  desired  to  tabulate  the  reduced 
data  then  enter  the  appropriate  record  numbers.  The 
reduced  data  will  be  tabulated. 
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TABLE  B-I 


VARIABLES  ADDED  TO  DATA  REDUCTION  PROGRAM 

A2  Lowest  Record  Number 

A 5 Counter 

A6  Decision  Variable 

A7  Highest  Record  Number 

B (Array)  Polynomial  Coefficients  From  TTR9 

F8  Decision  Variable 

Q (Array)  Values  of  Isentropic  Head  Coefficient 

T4  Polynomial  Approximation  of  Reynolds  Number 

Z Decision  Variable 


I 
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